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Abstract. - We study the effects of spatial constraints on the structural properties of networks 
embedded in one or two dimensional space. When nodes are embedded in space, they have a 
well defined Euclidean distance r between any pair. We assume that nodes at distance r have a 
link with probability p(r) ~ r~^ . We study the mean topological distance I and the clustering 
coefficient C of these networks and find that they both exhibit phase transitions for some critical 
value of the control parameter 5 depending on the dimensionality d of the embedding space. We 
have identified three regimes. When S < d, the networks are not affected at all by the spatial 
constraints. They are "small- worlds" I '^ logA^ with zero clustering at the thermodynamic limit. 
In the intermediate regime d < S < 2d, the networks are affected by the space and the distance 
increases and becomes a power of log A'^, and have non-zero clustering. When S > 2d the networks 
are "large" worlds I ~ A^^' '^ with high clustering. Our results indicate that spatial constrains have 
a significant impact on the network properties, a fact that should be taken into account when 
modeling complex networks. 



Introduction. — The interdisciplinary field of com- 
plex networks has recently received considerable atten- 
tion [1-7,9-12]. A network is a set of nodes, representing 
the elementary units of a system. These nodes are con- 
nected by edges (bonds) when an interaction between the 
nodes is present. In a social network, for example, the 
nodes denote individuals. An edge connects two nodes if 
there is a relationship between the two individuals. Net- 
works are particularly interesting as theoretical models of 
complex systems since they provide a straightforward way 
to describe all sorts of interactions between the elementary 
units of a system. 

Although there are infinitely many types of networks 
the attention of scientific community is actually focused 
in four types: Regular networks. Random graphs, "Small- 
world" networks and scale-free networks. 

Regular Networks consist of nodes that are identical to 
each other. The number of edges k emanating from each 
node (the "degree" of the node) is a constant. Some of the 
most common constructions used in computational and 
theoretical physics, such as the square lattice, the Bethe 
lattice or the complete graph (a graph where each node 
is connected to every other node) belong to this type of 
networks. 



Random graphs were introduced and studied by Erdos 
and Renyi [13]. In Erdos -Renyi networks the distribu- 
tion of the degrees of the nodes is Poissonian, and this 
fact facilitates the study of several important properties 
and allows analytical results. These networks have been 
extensively studied and used as some of the first models 
of social networks. 

"Small-world" networks were initially proposed by 
Watts and Strogatz as an improved model of social net- 
works [6,7]. These networks show simultaneously the im- 
portant feature of "clustering" (i.e. nodes that have com- 
mon neighbors have a higher probability of being them- 
selves connected) and the small- world phenomenon mean- 
ing that their topological diameter is slowly (logarithmi- 
cally) increasing with increased system size. 

In scale-free networks [1-3] finally, the degree distribu- 
tion follows a power law P{k) ~ k~'^ , with 7 typically 
between two and three. Scale-free networks have the ul- 
trasmall world property as the diameter of the network 
scales only logarithmically with the logarithm of the num- 
ber of nodes [10]. 

In all the above cases, the spatial arrangement of nodes 
is considered of no importance. It is assumed that spa- 
tial constraints can be neglected such that mean field 
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approaches become applicable. Many real networks are, 
however, embedded in 2D or 3D space. Examples are the 
Internet, airline networks, social networks, transportation 
and wireless communication networks [3, 14] etc. If these 
spatial constraints are important, mean field theory is no 
longer valid [15], and new approaches have to be devel- 
oped instead. Spatially constrained networks are a promis- 
ing type of complex networks whose importance has been 
only recently realized [16-19]. Up to now, most efforts 
deal with the simplest case, i.e. how to embed in the nor- 
mal Euclidean space complex networks with a given degree 
distribution, while allowing only closest neighbor connec- 
tions. If one allows additional long range connections then 
the resulting networks are more difhcult to study but much 
closer to the real complex networks. 

Here we study spatially constrained networks embed- 
ded in one and two dimensional space. We are interested 
in the simple case of networks that have a Poissonian de- 
gree distribution and where nodes are connected to each 



node i with probability P{r) 



„d-l^-S 



where c is deter- 



other with a probability p{r) 



i.e. depending on the 



Euclidean distance r between the nodes. The choice of a 
power law for the distance distribution is supported from 
recent findings on the topology of real complex networks, 
such as the Internet [20], human travels [21] and in social 
networks [22]. 

We propose a simple algorithm to generate such net- 
works and we study properties such as the topological dis- 
tance I between nodes and the clustering cocfhcient C. 
The exponent 6 has the role of a control parameter of 
the properties of the embedded networks and we find that 
there are three regimes of interest with qualitatively dif- 
ferent behavior. Networks with 6 < d, where d is the di- 
mensionality of the embedding space, behave more or less 
like random graphs, networks with S > 2d behave much 
like regular lattices while networks with S between these d 
and 2d show intermediate power-law behavior with critical 
exponents that depend on 6. 

Conceptually, our model is related to the "small-world" 
networks [6-8] and to long-range percolation [23-25]. In 
both the above cases, an underlying lattice structure is as- 
sumed and is "perturbed" by the addition of a fraction of 
long range connections. Thus, there is a competition be- 
tween short and long range interactions which determines 
the properties of the network. In our model, no such un- 
derlying structure exists. A perfectly ordered lattice is 
recovered only in the limit of infinite 6. There is, how- 
ever, still a "hidden" competition between short and long 
range, but now the magnitude of both "short" and "long" 
is controlled by the exponent S and the system size. 

Methods. — To construct the networks, we first ar- 
range the nodes in a one dimensional {d — I) linear chain 
or in a two dimensional (d = 2) regular square lattice. 
Thus, between any pair of nodes there is a well defined 
Euclidean distance. Next, we initially choose a node i ran- 
domly and assume that it is connected to kt other nodes. 
For each of these nodes, we first choose its distance r from 



mined from the normalization condition J. P{r)dr = 1, 
with L = N^/'^. Then we pick one of the Nr sites that 
are at distance r within the underlying lattice. We repeat 
this process for all nodes i in the underlying lattice and 
then remove multiple connections. The nodes of the final 
network we obtain this way do not have a fixed degree. 
Their mean degree k depends on both kt and 5. In other 
words, for obtaining a certain k value for fixed 5, we have 
to choose an appropriate value for kt . When the networks 
are constructed we measure their actual degree probability 
distribution p{k) to verify that all networks have at least 
approximately the same mean degree k and the actual dis- 
tance distribution p(r) i.e. the fraction of connected nodes 
that are at distance r from each other, which is affected 
by the process of the removal of duplicate edges. In what 
follows we have succeeded in obtaining fc « 4 (to be pre- 
cise 3.3 < fc < 4) using kt = 2 when i5 < 2.25, kt = i when 
2.25 < (5 < 3 and fcj = 5 for (5 = 4 for the one dimensional 
case. For the two dimensional case kt — 2 produced net- 
works with 3.6 < fc < 4 for all 5 < 4.5. Thus, no particular 
"fine-tuning" was necessary. 

Quantities of interest. — We are interested in the 
mean topological distance / between the nodes of the net- 
work and the mean clustering coefficient C. The topolog- 
ical distance between two nodes a and b is the minimum 
number of links that a walker has to cross in order to ar- 
rive from a to b. The (local) clustering coefficient of a 
node I in a network is defined as: 



C{h) = 



U 



[h{h - l)/2] 



(1) 



where ti is the number of triangles (loops of length 3) 
attached to this node divided by the maximum possible 
number of such loops. For calculating the mean topolog- 
ical distance I and the mean clustering coefficient C we 
average over all nodes in the network (typically 10"'') and 
over all network (typically 10^ ) realizations. 

Results. — First, we show that the generated net- 
works have the desired degree and distance distributions. 
Figure 1 shows the degree distribution P{k) and the dis- 
tance distribution P{r) for the generated one dimensional 
(Fig l(a,c)) and two dimensional (Fig l(b,d)) networks, for 
three characteristic 5 values each, S ~ 0.5, 1.5, 2.5 {d = 1) 
and S = 2,3.5,5.5 {d = 2). Figures la,b confirm that the 
networks have a narrow degree distribution around A: « 4. 
Figures lc,d confirm that P{r) has the desired distance 
distribution P{r) - r-^'^-'^+'^'> for both d == 1 and rf = 2. 

Next, we study the mean topological distance I as a 
function of the network size N. In the limiting case (5 ^- 
the spatial constraints are not significant and we expect 
that the networks will belong to the universality class of 
Erdos -Renyi random graphs, where I scales as I ~ log TV. 
In the opposite case S —>■ oo, the spatial constraints will 
result in an ordered lattice-like structure that belongs to 
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10 10 10' 10' 10 10 10 10" 10" 
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Fig. 1: (a) Probability P{k) tliat a node has degree k for d — 1 
networks with A^ = 10^, fc fa 4 and for S = 0.5,1.5 and 2.5. 
(b) Probability P{k) for d = 2 networks with Af = 5 x lO", 
fc « 4 and for S — 2, 3.5, 5.5. (c) Probability P{r) that a node 
has a connection at distance r for the same d = 1 networks 
as in (a). The points are logarithmically binned simulation 
data, (d) Probability P{r) that a node has a connection at 
distance r for the same d = 2 networks as in (b). The points are 
logarithmically binned simulation data. In (c,d) the straight 
lines have slopes equal to — (5 — d + 1). 




N(size of the largest cluster) 

Fig. 2: (a)Mean topological distance I as a function of the net- 
works size A'^ for d = 1 networks with A; « 4 and for S = 0, 0.5 
and 1 i.e. in Regime 1 (LinearLog plot). (b)Mean topologi- 
cal distance I as a function of the networks size N for d — 2 
networks with A; « 4 and for 5 — 0,1 and 2 i.e. in Regime 
1 (LinearLog plot). The straight line is guide to the eye. (c) 
Mean topological distance / as a function of the networks size 
TV for d = 1 networks with fc ~ 4 and for S — 2, 2.50 and 4 
i.e. in Regime 3(LogLog plot).(d) Mean topological distance / 
as a function of the networks size A*' for d — 2 networks with 
fc ~ 4 and for S — 4,5, 5.5 i.e. in Regime 3 (LogLog plot). The 
straight line has slope 0.5. 



the universality class of regular lattices, where the mean 
topological distance I scales as N^^'^ ( d is the dimension- 
ality of the embedding space). We are interested in how 
the dependence of the topological distance I on N changes 
when S changes from zero to infinity. 

Our results (Figs. 2 and 3) indicate that in both d = 1 
and d = 2, there are three regimes separated by ^ = d and 
S = 2d, where the behavior of the networks is qualitatively 
different. In the first regime < 6 < d (see Fig. 2a for 
d — 1 and Fig. 2b for d = 2), I scales as logiV independent 
of 6. In the third regime, 6 > 2d (see Fig. 2c for d = 1 and 
Fig. 2d for d = 2) the result suggests that I asymptotically 
scales as iV^/'* as for regular lattices. This dependence is 
well observed in d = 1 (Fig. 2c) for (5 = 4 and in d = 2 for 
S ~ 5.5. For 5 closer to 2d, this scaling is only recovered 
for sufficiently large networks. For S = 2d, I scales as N'', 
with 7 « 0.35 for d = 1 and 7 w 0.3 for d = 2 (see the 
lowest curves in Figs. 2c and d, respectively). 

Figure 3(a,b) shows the mean topological distance I as a 
function of N in the intermediate reginre d < S < 2d. The 
plots suggest that for both d = 1 and d = 2,1 scales as I ~ 
(ln(7V))". Accordingly, we can sunrmarize the dependence 
of Z on A^ as follows: 



log(iV) 

(^og{N)Y 



S <d 
d<S <2d 
S>2d 



(2) 



When S = 2d the numerical results are quite well described 
by a power law A^'*' where 7 w 0.35 for d = 1 and 7 « 
0.3 for d = 2. It is, however, probable that there are 
logarithmic corrections to this scaling form. The variation 
of the exponent a with d is shown in Fig. 3c for d = 1 and 
Fig. 3d for d = 2. The results suggest that approximately 



a 




(3) 



Equation (3) provides a good fit to the results in Figs. 
3(c,d), and although is not a rigorous analytic result it may 
be useful for practical purposes ^ Note that Eq. (3) ensures 
a smooth transition between the first and second regime. 
To match with the power law behavior when S > 2d, a 
diverges when approaching S = 2d fronr below. In d = 2, 
also a weaker singularity could fit the data. 



^Of course, in Eq. (2) in the range d < 5 < 2d st more complex 
scaling form cannot be excluded. 



p-3 



Kosmas Kosmidisl, Shlomo Havlin 2Armin Bundel 




1.5 

Exponent § Exponent 8 

Fig. 3: (a) Mean topological distance I as a function of the net- 
works size A'^ for d = 1 networks with fc ~ 4 and for S — 1.2, 1.5 
and 1.8 (bottom to top) i.e. in Regime 2 (LinearLog plot). The 
lines are fits to the simulation data of the scaling assumption 
(Eq. 2).(b) Mean topological distance I as a function of the net- 
works size A^ for d = 2 networks , A; « 4 and for S — 2.5, 3, 3.5 
(bottom to top) i.e. in Regime 2 (LinearLog plot). Points 
are simulation data and lines are fits of Eq. 2.(c) Exponent 
Q as a function of S for d = 1 networks in the second regime 
(1 < (5 < 2). (d) Exponent a as a function of S for d — 2 
networks in the second regime (2 < 5 < 4). Points are best fit 
estimates to the simulation data. Lines are plots of Eq.(3) 



We like to note that these transition points arise natu- 
rally in the scaling of the mean distance f and the maxi- 
mum distance rmax with the system size L. By definition 
r = J-, rp{r)dr, while r^nax is obtained from the condition 
L'^ J p{r)dr w 1. It is easy to verify that 



L S <d 

where ci does not depend on A^ and 



L 



0<S <2d 
6>2d 



(4) 



(5) 



Accordingly, the relative mean distance f/L is finite for 
S < d and for S > d in the limit of large L. Similarly, 
the relative maximum distance rmax/ L is finite for 5 <2d 
and for 5 > 2d for L going to infinity. The above rela- 
tions can qualitatively elucidate the existence of different 
regimes. For simplicity consider the case d — \ and as- 
sume one attempts to go from the one end of the "chain" 
to the other. When the mean length f of each link is 
comparable to the lattice size L then he can succeed with 
only one step on average. The lattice is a small world 
(regime 1). When the the mean length f is less than L 
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is of the order of L, then the walker will have to 



Fig. 4: (a) Clustering coefficient C as a function of 1/A^ for 5 = 
0,0.5, 1, 1.2, 15, 1.8, 2.5 and d = 1. (b) Clustering coeflicient C 
as a function of IjN for 5 = 1, 1.5, 2, 2.5, 3, 3.5, 5 and d = 2. 
Straight lines are best fits with slope 1. 



do more steps but he will eventually arrive at a site that 
has a long edge r^nax that will allow him to arrive at the 
end of the lattice. The lattice is not a small world as in 
the first case but also not a regular lattice since there are 
links that decrease the topological distance considerably 
(regime 2). When both f and r^ax are less than L, there 
are no shortcuts. Renormalizing the graph will result in a 
regular linear chain (regime 3). 

Finally, we study how the clustering coefficient C de- 
pends on 6 and A^ (Fig. 4). We are particularly interested 
in the behavior at the transition points 5 — d and 5 = 2d. 
We expect that for 5 — > 0, C ~ 1/A^ similar to the case of 
Erdos -Renyi random graphs. When S ^ oo we expect the 
(high) clustering of the lattices. We find that the critical 
point ai S = d separates a phase of zero clustering when 
5 < d from a phase of finite clustering when S > d. The 
second transition point at d = 2d does not influence C re- 
markably. Figures 4(a,b) show C as a function of 1/A^ for 
d = 1 and d ~ 2 respectively. The clustering coefficient C 
decreases as 1/A^ for S < d in both one and two dimen- 
sions, while for (5 > d the slope of the curves in the double 
logarithmic plot becomes zero, confirming the existence of 
a regime with non zero clustering. It is remarkable that as 
soon as 6 becomes larger than d, the spatial effects become 
very pronounced. The topology not only affects the mean 
topological distance but also increases simultaneously the 
clustering coefficient. 

Conclusions. — We have examined the mean topo- 
logical distance I and the clustering coefficient C of net- 
works that are embedded in one or two dimensional space. 
Both exhibit phase transitions for some critical value of the 
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control parameter 6. We have identified three regimes. 
When S < d, where d is the dimensionality of the em- 
bedding space, the networks are "small-worlds" with zero 
clustering at the thermodynamic limit. In the interme- 
diate regime d < S < 2d, the topology changes. The 
topological distance becomes a power of the logarithm of 
the system size and the networks have non-zero clustering. 
When S > 2d the networks do not have the small world 
property and have high clustering. Our results indicate 
that spatial constraints have an important influence on 
the network properties and this fact should be taken into 
account when modeling real complex networks. A natural 
extension of this work will be to consider fractal or mul- 
tifractal substrates instead of regular lattices [14, 26] and 
study their influense on the topological properties of the 
embeded networks. 
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